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Abstract 

Background: Group 1 introns (ribozymes) are among the most ancient and have the broadest phylogenetic 
distribution among the known self-splicing ribozymes. Fungi are known to be rich in rDNA group 1 introns. In the 
present study, five sequences of the 28S ribosomal RNA gene (rDNA) regions of pathogenic dematiaceous 
Phiolophoro verrucosa were analyzed using PCR by site-specific primers and were found to have three insertions, 
termed intron-F, G and H, at three positions of the gene. We investigated the distribution of group 1 introns in this 
fungus by surveying 34 strains of P. verrucosa and seven strains of Phialophora americana as the allied species. 

Results: Intron-F's (inserted at L798 position) were found in 88% of P. verrucosa strains, while intron-G's (inserted at 
LI 921 ) at 12% and intron-H's (inserted at L2563) at 18%. There was some correlation between intron distribution 
and geographic location. In addition, we confirmed that the three kinds of introns are group 1 introns from results 
of BLAST search, alignment analysis and Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR). Prediction of 
secondary structures and phylogenetic analysis of intron sequences identified introns-F and G as belonging to 
subgroup IC1. In addition, intron-H was identified as IE. 

Conclusion: The three intron insertions and their insertion position in the 28S rDNA allowed the characterization 
of the clinical and environmental isolates of P. verrucosa and P. americana into five genotypes. All subgroups of 
introns-F and G and intron-H were characterized and observed for the first time in both species. 
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Background 

The type species Phialophora verrucosa was described by 
Medlar in 1915 [1] when he isolated the fungus from a 
human skin disease. The species is ubiquitous and cos- 
mopolitan, and are important plant saprobes as well as 
human pathogens. Identification is based on conidial 
ontogeny and molecular systematics. Few studies invol- 
ving molecular genotyping techniques have been 
reported for P. verrucosa. A study analyzed restriction 
fragment length polymorphisms (RFLP) of mitochondrial 
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DNA to determine genetic variations and phylogenetic 
relationships among P. verrucosa strains [2] . 

Different molecular typing tools, such as random 
amplification of polymorphic DNA (RAPD), RFLP, 
pulsed-field gel electrophoresis (PFGE), multilocus 
enzyme electrophoresis (MLEE) and multilocus 
sequence typing (MLST), have been developed to pro- 
vide a better understanding of the molecular epidemiol- 
ogy of fungal pathogens, e.g., Candida albicans [3-5] 
and Aspergillus fumigatus [6,7] and medically important 
filamentous fungi [8]. However, although the majority of 
the reported group 1 intron sequences have been found 
in a wide range of fungi (Comparative RNA Web 
[CRW] site: http://www.rna.ccbb.utexas.edu/l9], few stu- 
dies about sequence and structure variation, distribution 
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and phylogenetic relationships of introns from a single 
species have been performed in detail. We focused on 
group 1 introns within 28S rDNA from P. verrucosa to 
evaluate the prevalence of intron polymorphism at the 
strain level. 

As the first step to determine intron sequence diver- 
gence, sequences of 28S rDNA of five representative 
strains of P. verrucosa were analyzed to find insertions. 
Based on these five sequences, site-specific primers were 
designed for use in PCR to detect insertions on other P. 
verrucosa and P. americana strains studied, in order to 
investigate incidence and distribution of insertions. 
Moreover, to characterize the insertions, we analyzed 
the phylogeny of the introns found in this study and 
predicted their secondary structures. 

Results 

Nucleotide structure of P. verrucosa 28S sequences and 
the characterization intronic insertion 

Since the sequence information of the 28S region was 
not available in public databases, we sequenced the 
28S, including ITS, regions of five representative 
strains (Table 1) of P. verrucosa. Alignment of the five 
P. verrucosa rDNA sequences revealed the nucleotide 
sizes of 18S, ITS and 28S and the location of intron-F 
and G insertions (see Table 2 and Additional file 1). It 
was found that the sequences were composed of 29 
nucleotides of partial 18S, 534-535 nucleotides of ITS 
and from 3349 to 4133 nucleotides of 28S regions. 
These genomic sequences were deposited at DDBJ and 
accession numbers are listed in Table 2. Twenty-five 
and two nucleotide substitutions were found within 
the ITS and D1D2 regions, respectively. The two bp 
substitutions in D1D2 were located at 1036 and 1042 
nucleotide positions in the 28S region. These poly- 
morphisms were confirmed in the five strains of P. ver- 
rucosa strains except for the two insertions in intron-F 
and G at 924 and 2239 positions of the 28S, respec- 
tively. None of the insertions were found in the Yao 
strain. BLAST comparisons showed with relatively high 
homology values that the sequences of both insertions 
were individually homogeneous to group 1 introns in 
the database. The positions of insertion for intron-F 
and G correspond to nucleotide positions 798 and 
1921 of Escherichia coli 23S nucleotide sequence 
accession number J01695 [10]. In addition, although 
these five strains of P. verrucosa did not possess 
intron-H, this region of PV28 strain was amplified 
with the site-specific primer pair (in Table 3) designed 
from the sequences obtained from another experiment 
(data not shown) and a 403-nucleotide insert repre- 
senting intron-H was sequenced. The insertion posi- 
tions of the intron-H correspond to nucleotide 
positions 2563 of E. coli 23S nucleotide sequence. 



Survey of insertions of P. verrucosa and P. americana 

We amplified intron insertion regions using site-specific 
primer pairs we have designed for intron-F (inF-F and 
inF-R), intron-G (inG-F and inG-R) and intron-H 
(L2563F and L2563R), within the 28S region (Table 3). 
These primer pairs were used to screen and detect PCR 
amplicons for insertion regions within 34 P. verrucosa 
and seven P. americana strains. Amplicons were eluted 
in agarose gel to gain information regarding the intron 
insertions. No-insertion amplicons for intron-F and 
intron-G primers were in the size 142 andl85 bps, 
respectively. When insertions were present, intron-F pri- 
mer pair yielded amplicons in the size range from 531 
to 533 bps, and intron-Gs in the size 575 or 578 bps. 
Moreover, amplicons of about 643 bps for intron-Hs 
were also eluted. It was revealed that there were 30 
intron-F's, four intron-G's and six intron-H's within P. 
verrucosa and only two intron-Fs within P. americana 
as shown in Table 1. There was some correlation 
between intron distribution of P. verrucosa and geo- 
graphic location, i.e., intron-Fs were found to have pre- 
valence of 88% in P. verrucosa and intron-Hs were 
found specifically in the South American Continent. No 
introns were found except for two intron-Fs in P. ameri- 
cana. In addition, the agarose gel profiles allowed us to 
characterize genotypes and distribution frequencies of 
insertions from P. verrucosa including no-insertion as 
shown in Table 1. It was found that occurrence of geno- 
types F, FG, FH, FGH and N were at 64, 6, 12, 6 and 
12%, respectively. 

Characterization of the P. verrucosa intronic insertion 

RT-PCR was carried out to identify the property of 
these insertions, namely, whether they are introns or 
unusual extensions incorporated into mature rRNA. 
Four representative strains were selected among the 41 
strains surveyed. And it was found that two strains (PV1 
and PV3) had two introns individually, while the other 
two strains (PV2 and PV41) had only one intron as 
shown in Figure 1. Insertions of strain PV1 and PV3 
were eluted at 142 bps on lane 2 and 3 with intron-F 
primer pair, and 185 bps on lane 4 and 5 with intron-G 
primer pair, respectively. PV2 and PV41 exhibited 142 
bps amplicons with intron-F primer pair as shown on 
lane 15 and 16, respectively. An intron-lacking Yao 
strain gave 142 and 192 bps amplicons with intron-F 
and G primer pairs on lane 10 and 11, respectively. The 
other lanes; namely, 6, 7, 8, 9, 13 and 14 show PCR pro- 
ducts of genomic DNA as templates and lane 12 is 
negative control. These shows that all the insertions 
were excised after cDNA was transcribed, indicating 
that they were, as predicted, actively spliced introns. 
These results point to the possibility that these inser- 
tions are group 1 introns. 
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Table 1 Thirty-four strains of P. verrucosa and seven strains of P. americana examined in this study 


Strain 


Isolation source 


Locality 


Insertion 


Genotype 


Sample ID 








Intron-F Intron-G Intron-H 






IFM 4928 


Human 


Japan 


+ + - 


FG 


PV1 


IFM 5089 


Human 


Japan 


+ - - 


F 


PV2 


IFM S1934 

II 1 VI J 1 yjT 


Hi imsn 
i i u i i la i i 






N 


Yao 


IFM 41 77Q 

II IVI t 1 / / 7 


Hi iman 
i i u 1 1 id i i 


China 


+ + 


FG 


PV3 


IFM 41 7SD 

II IVI T I / JU 


Bark 


China 


+ 


F 


TH9 


IFM41 71 D 

II IVrr 1 / 1 U 




China 


+ - na 


F (na) 
r {\ \d) 




IFM 41 7R0 

II IVI "T I / ou 


Hi iman 
i i u 1 1 id i i 


China 




F 




IFM 41 7?1 

II IVI tl / Z. 1 


Rnttpn \A/nnH 


China 


_i_ _ _ 


F 




IFM 41 774 

II IVI T" I / 


Rnttpn \A/nnH 
RULLCI 1 vvuuu 


China 


+ 


F 




IFM 41 7?S 

II IVI T" I / Z..J 




China 


+ 


F 




IFM 41739 


Bark 


China 


+ - - 


F 




IFM 41740 


Bark 


China 


+ 


F 




IFM 41746 


Straw 


China 


+ 


F 




IFM 41749 


Bark 


China 


+ - - 


F 




IFM 41752 


Bark 


China 


+ - - 


F 




IFM 41764 


Soil 


China 


+ 


F 




IFM 41755 


Unknown 


China 


+ 


F 




IFM 41765 


Soil 


China 


+ 


F 




IFM 41778 


Straw 


China 


+ 


F 




IFM 41872 


Soil 


Colombia 




N 




IFM 41871 


Soil 


Colombia 


+ - - 


F 


PV41 


IFM 41879 


Soil 


Colombia 


+ - + 


FH 




IFM 41881 


Soil 


Colombia 


+ 


F 




IFM 41885 


Soil 


Colombia 


+ 


F 




IFM 41884 


Cactus 


Venezuela 




N 




IFM 41888 


Rotten wood 


Venezuela 




N 




IFM 41887 


Rotten wood 


Venezuela 


+ - + 


FH 


PV28 


IFM 41886 


Cactus 


Venezuela 


+ - + 


FH 


TH31 


IFM 41892 


Soil 


Venezuela 


+ - + 


FH 


TH35 


IFM 41883 


Soil 


Venezuela 


+ 


F 




IFM 41893 


Soil 


Venezuela 


+ 


F 




IFM 41897 


Soil 


Brazil 


+ + + 


FGH 


PV33 


IFM 41 RQR 
Irivi H- 1 oyo 


jUN 


Rravil 
Dl dZM 


+ + + 


Fr,H 

run 


P\/34 
r V dh 


IFM 41899 


Soil 


Brazil 


+ 


F 




*CBS 273.37 


Human 


Brazil 




N 




*CBS 400.67 


Soil 


Brazil 




N 




*CBS 281.35 


Human 


USA 




N 




*CBS 220.97 


Linden tree 


USA 




N 




*CBS 840.69 


Decaying timber 


Finland 




N 




*CBS 221.97 


Unknown 


Uruguay 


+ 


F 




*CBS 223.97 


Human 


USA 


+ 


F 





*: P. americana, +: with insertion, -: no insertion, na: not analized. 



Moreover, we analyzed sequences of the spliced 
introns to confirm the boundaries of exon and intron 
sequences. The last nucleotide of the upstream exon 
was confirmed to be a T (U in RNA) and the last 
nucleotide of the intron was a G, consistent with group 
1 introns [11,12]. 



Phylogenetic relationships of introns F and G of P. 
verrucosa 

Sequences of intron-F and G of ten P. verrucosa strains 
were sequenced and it was found that DNA sequence 
polymorphisms exist among the two introns, i.e., the 
intron-Fs ranged in the size from 389 to 391 bps and 
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Table 2 List of ITS, 28S rDIMA and intron sequences of P. verrucosa 



Sample ID or entry name 






Length (bp) 


Splice positions 


Accession number 


ITS 


28S 


Intron-F Intron-G Intron-H 


position 3 


position' 3 




PV1 


535 


4130 








AB550775 


PV2 


535 


3922 








AB550776 


PV3 


535 


4133 








AB550777 


PV41 


534 


3922 








AB550778 


Yao 


535 


3349 








AB550779 


F-PV1 






391 


924 


798 




F-PV2 






391 


924 


798 




F-PV3 






391 


924 


798 




F-PV41 






391 


924 


798 




G-PV1 






390 


2239 


1921 




G-PV3 






393 


2239 


1921 




F-TH9 






389 


924 


798 


AB550780 


F-PV28 






389 


924 


798 


AB550781 


F-TH31 






389 


924 


798 


AB550782 


F-TH35 






389 


924 


798 


AB550783 


F-PV33 






390 


924 


798 


AB550784 


F-PV34 






390 


924 


798 


AB550785 


G-PV33 






389 


2239 


1921 


AB550786 


G-PV34 






389 


2239 


1921 


AB550787 


H-PV28 






403 


2905 


2563 


AB611046 



a Position means relative to the 28S rRNA of P. verrucosa Yao strain and b position means relative to 23S rRNA of E. coli J01965. 



the four intron-Gs from 389 to 393 bps shown in Table 
2. There were 24 nucleotide substitutions and two dele- 
tions/insertions (TH9 strain) within intron-F. There 
were five nucleotide substitutions among intron-Gs 
from PV1, PV33 and PV34, unlike 36 substitutions 
between PV1 and PV3. In addition, Blast search analyses 
and alignment lead us to believe that intron-Fs and Gs 
from 14 introns belong to subgroup IC1 of group 1 
intron. 

Fourteen introns from 12 representative strains of P. 
verrucosa including Tetrahymena thermophila as out- 
group were aligned and used for phylogenetic analyses. 
Neighbor-joining (NJ) and Maximum Parsimony (MP) 
trees based on the alignment of these intron sequences 
are shown in Figure 2. The data set consisted of 466 
characters, of which 156 were removed from the MP 
analysis due to ambiguous alignment. Of the remaining 
310 characters, 201 were variable and 129 were phylo- 
genetically informative for parsimony analysis. Three 
major distinct and well-supported clades that had 
homologous topology were obtained from both phyloge- 
netic analysis methods showing that all the introns ana- 
lyzed were undergoing a similar rate of evolution. The 
first clade [I] (87% BS support in NJ, 81% in MP) con- 
sisted of six strains having intron-F including 3 clinical 
isolates, the second clade [II] (57% BS in NJ and 77% in 
MP) consisted of 4 strains having intron-F, and the 
third clade [III] (100% BS in both trees) consisted of 



four G introns. All the introns clustered in clades [I] 
and [II] are inserted at the same position L798 those in 
clade [III] at the same position L1921. Introns inserted 
at the same positions belong to the same clusters and 
are considered to be the same subgroups. We observed 
that two Venezuelan strains in clade [II] (TH31 and 
TH35) included the remnants of homing endonuclease 
(HE) gene sequences, which were found to invade into 
IC1 introns of the 18S regions in another study (unpub- 
lished data). 

Alignment and phylogenetic analysis of the core regions 
of the group IC1 introns 

Alignment of the core regions consisting of highly con- 
served sequences of the elements of P, Q, R and S and 
the pairing segment P3 and the nucleotide sequences, in 
particular, the last two nucleotides GC of the Q element 
and the first and second GU nucleotides of the R ele- 
ment [12] (Additional file 2) showed that the introns 
belong to group IC1. All core region sequences of 
intron-Fs were found to be identical. Two sequences of 
core regions termed as intron-G (PV3) and intron-G 
(PV1, PV33, PV34) were obtained and added to the NJ 
analysis in Figure 3. The NJ tree was constructed based 
on the alignment of these core regions consisting of 
three representative sequences of P. verrucosa and IC1 
of 21 taxa drawn from database using IE intron from 
Neoscytalidium dimidiatum as out-group. The 
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Table 3 Primers used for the amplification and sequencing of P. verrucosa 


Primer 


Sequence (5'-3') 


5' position* 


Source 


5' position including ITS 


ITS1 


TCCGTAGGTGAArfTGCGG 


-563 


Whitp Tl pt al T4R1 

VVIIIIC I J, CI CM. L'^J 


1 


ITS 3 


G,G ATGG, ATG, A A G, A A GG,G A G,G 


-309 


Whitp Tl Pt al 

Willie I J, cL al. LiOJ 


255 


NL1 


r;r AT ATfAATAAGCGGAGGAAA 


39 


O'Donnpll K T491 

w L^Vyl II ICTII l\ |_ i-^J 


603 


3PV26 


CCGTCTTGAAACACGGACC 


633 


Thic wnrk 

I 1 llo VVVJI l\ 


1 1 97 


inFG-F 


rrc^ a a ag atc^tc a agtatg,gg 


795 


Thic wnrk 

1 1 llo VVVJI l\ 


1 359 


inF-F 


A G GYGC A A ATfG AJGGJG A A 

AA^.vJ 1 vJ^.AAAAAA 1 LUn 1 ^.vJ 1 *^.AAAA 


868 


Thic work 

1 1 II J VVKJI l\ 


1432 


inF-R 


GAAGiGiGGTGTAATGATTGGiGT 

^_AAAAvJvJ^_^_ 1 ^_ 1 AAAA 1 ^_AA 1 1 1 


1 009 


Thic work 

1 1 II J VVKJI l\ 


1 573 


8PV26 


r, a a r rTTTr G G G A GTTG A G 

vJAAAA^_^_ 1 1 1 ^_^_^_^_AA^_ 1 1 ^_AAvJ 


1487 


Thic work 

1 1 II J VVKJI l\ 


2051 


1 1 PV26 


AAnrrATA^Q^AAGTTrnrr 

AAAAvJ^_^_AA 1 / \\3\3\3/\r\\3 1 1 V_V_V_J 1 


1 525 


Thic wnrk 

1 1 llo VVVJI l\ 


2089 


9PV26 


GTCGTACTC AT A AG 'CGCAG 

\J 1 ^.VJ | AAv_ 1 V^AA 1 AAAA^_^_vJ^_AAvJ 


181 8 


Thic \hirwV 

1 1 llo VVVJI l\ 


2382 


CA-INT-L 


ATAAGGGAAGTCGGCAAAATAGATrCGTAA 

n 1 AAAAVJvjvJAAAAVJ | \_\J\3^./\/\/\r\ 1 AAVJAA 1 1 AAAA 


1 881 


MrGi illni mh M 1 pt al TSD1 
i viv.v_u i iv^uy i i ivij, cl ai. L-'^J 


2445 


2PV26 


TrrrfiAAGTTArnGATrTA 

1 V_V_V.vJAAAAvJ 1 1 AAV_vJVjAA 1 V_ 1 AA 


1 91 8 


Thic work 

1 1 llj VV(J\ l\ 


2482 


1 6PV26 


GGG AAGGi TTACnACnCCA ATC~ 

^_V_V_AAAAV_V_V_ 1 1 AAvJAAvJV_^_AAAA 1 V. 


1 942 


Thic work 

1 1 llj VV(J\ l\ 


2506 


1 0PV26 


GGGJAGGA GG\JGT A A GJTG 

^_V_vJ 1 AAV_V_AAvJ 1 1 L 1 AAAAvJ 1 1 vJ 


2089 


Thic \hirwV 

1 1 1 lO VVVJI l\ 


2653 


inG-F 


GATGGCCAGAAAGTGGYGYTG 

vJAA 1 vJvJ^_V_AAvJAAAAAAvJ 1 vJvJ 1 VJ I I VJ 


21 30 


Thic wnrk 

1 1 1 lO VVVJI l\ 


2694 


inG-R 


t a (irzrz A G A GJG,G,G, A A TGTG GT 

1 AAV3VjVJAAV_AAVJ 1 VJVJVJAAAA 1 V_ 1 V_VJ 1 


2314 


Thic wnrk 

1 1 1 lO VVVJI l\ 


2878 


i.UJ 1 1 N I J 


GTAGrGAAArrAfAGrrAAG 

V_ 1 AAVJV_VJAAAAAAV_V_AAV_AAVJV_V_AAAAVJ 


2323 


Thic work 

1 1 1 1 J VVVyl l\ 


2887 


(~A-I NT-R 


rrTTGGrTGTGGTTTrGrTAGATAGTAGAT 

V_V_ 1 1 VJVJV_ 1 VJ 1 VJVJ 1 1 1 V_VJV_ 1 AAVJAA 1 AAVJ 1 AAVJAA 1 


2343 


MrGi illni mh M 1 pt al TSD1 
i viv.v_u 1 1 v^uy i i ivu, ci cm. L-'^'J 


2907 


inFG-R 


GrTrTrrrArrTATTrTAGArr 

OV_ 1 V_ 1 V_V_V_AAV_V_ 1 AA 1 1 V — 1 AAV_AAV_V_ 


2445 


Thic xA/nrlc' 

1 1 ID VVUI l\ 


3009 


12PV26 


TGGTAmCACCGGCGATO 


2464 


This work 


3028 


IGS1-M-R 


ctgccacaagccag™tc 


2805 


This work 


3369 


IGS3-M 


OTCGATGTCGGCTOTCCT 


2838 


This work 


3402 


IGS-L 


TAGTACGAGAGGAACCGT 


2991 


Williamson ECM et al. [51] 


3555 


3IGS-PV 


TCTAAGTCAGAATCCGTGCCG 


3090 


This work 


3654 


5IGS1-PV 


ACGAGCTACTGAGCGTAAG 


3318 


This work 


3882 


6IGS-PV 


GACCACAGTCAGGCTOCG 


3349 


This work 


3913 


L2563 F 


CACAGGGATAACTGGOTGTGG 


2781 


This work 


3345 


L2563R 


ATCTGAATCAACGGTOCTCTCG 


3018 


This work 


3582 



* The 5' position is relative to the 28S rDNA sequence of the P. verrucosa Yao strain. 



phylogeny of intron-F and G formed separate clades as 
shown in Figure 3, and indicated that both introns were 
likely acquired independently. Indeed, all intron-Fs were 
found to be closely related to Myriosclerotinia ciborium 
and Sclerotinia tetraspora introns which are located at 
L798. Two sequences of intron-G located at L1921 were 
grouped together with 85% BS value and found to be on 
the neighboring clade with Cordyceps prolifica intron 
located at L1921. The phylogenetic tree suggests that 
both introns may be inserted prior to the divergence of 
the species formerly belonging to clade [IV] and [V]. 
Collectively, this tree displays that all introns of P. ver- 
rucosa generated by the core regions are members of 
subgroup ICls. 

Modeling of the P. verrucosa insertions revealed that they 
were group IC1 introns 

The predicted secondary structure of the intron-F and G 
were constructed as follows. The conserved P, Q, R and S 
regions of intron-F (L798) from PV1 were initially aligned 
with the same regions from other taxa, and then regions 



of PI through P10 were constructed and added on the 
basis of the secondary structure model as shown in Figure 
4[A][13]. The P10 region was formed using the internal 
guide sequence (IGS). The catalytic core was defined by a 
set of structurally conserved elements, including elements 
P3 to P8. A G-C pair within P7, i.e. G391-C277 of intron- 
F was assumed to be G-binding positions [14]. Extended 
P5 and P9 stems were displayed in the putative structure 
of intron-F from PV1. Nine intron-Fs from nine strains 
(PV2, 3, 28, 33, 34 and 41 and TH9, 31 and 35) of P. verru- 
cosa were predicted to be the same structures as the puta- 
tive structure of intron-F derived from PV1 drawn in 
Figure 4[A], alternatively, shown in Additional file 3. 
These nucleotide variations among intron-F were observed 
mainly in the loop and at four positions where one nucleo- 
tide of P5a, two of P5.1a and one of P5.2 stem were posi- 
tioned. The base pairs GU and CG within P6 were formed 
in the core region of intron-F [12]. The nucleotides A71, 
A72, U73 were located in segments J3/4 of PV1 intron-F 
[15-18]. These predictions of secondary structure revealed 
that all intron-Fs were IC1 group 1 introns. 
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1,000bp 




Figure 1 Amplification pattern by RT-PCR with the site-specific 
primer pairs for intron-F and G. PCR products of from cDNA 
amplified with the primers inF-F and inF-R are eluted in lanes 2, 3, 
15 and 16, and with primers inG-F and inG-R in lanes 4 and 5. PCR 
products from genomic DNA amplified with primer pair for intron-F 
are eluted in lanes 6, 7, 10, 13 and 14, and with primer pair for 
intron-G in lanes 8, 9 and 1 1. Lane 12 is the negative control. 



The structure of intron-G (L1921) from PV1 was 
drawn just as was done for intron-Fs (Figure 4[B]). A G- 
C pair within P7, i.e. G390-C360, was assumed to be the 
G-binding positions. The GU-CG pair of P6 and the 
AAU in J3/4 was the same as in the intron-F core 
region of PV1. This putative intron-G exhibited 
expanded regions of PI and P5. The three intron-Gs of 
PV1, PV33 and PV34 were found to be similar among 
the three strains. Different features were found in PV3 
as shown in Figure 4[C] wherein the sequence of PV3 
differed in PI region among four trains; namely, short 
stems in Plb and Pic and small bulge loops of LI and 
Lla (Additional file 4). Moreover, PV3 added P2.0 and 
P8c, although the other intron-Gs did not. Prediction 
structures in the remaining two introns of PV33 and 
PV34 are not shown. Nevertheless, all subgroups of 
intron-G were also identified as IC1, based on 
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Figure 2 Phylogenetic relationship of intron-F and G within 

28S of P. verrucosa. The trees were generated using MP (A) and 
NJ (B). One of three equally MP trees (tree length = 353, 
consistency index (CI) = 0.9575, homoplasy index (HI) = 0.0425, CI 
excluding uninformative characters = 0.9268, HI uninformative 
characters = 0.0732, retention index = 0.9679, rescaled consistency 
index = 0.9268). * indicates a clinical isolate of P. verrucosa. 
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Figure 3 Phylogenetic tree of IC1 intron based on elements P, 
Q, R, S and a segment of P3. Numerals at each node are 
bootstrap probabilities from NJ analysis. Insertion positions are given 
after the sample ID or accession number. * indicates the insertion 
position relative to the 18S rDNA of the 5. cerevisioe sequence. 



comparison of tertiary structures across segments P3-7 
of the four strains. 

In conclusion, we have identified that the ten intron- 
Fs and four intron-Gs of P. verrucosa belong to IC1 
group 1 introns. 

Characterization of intron-H 

Loss of P5abcd domain in derived S788 introns was cor- 
related with inability to self-splice in vitro in a previous 
report [19]. Accordingly, we have not confirmed inser- 
tion positions of intron-H by RT-PCR. However, we 
examined PV-28 strain as the representative strain of 
intron-H by analyzing the sequence alignment of the 
core region of subgroup IE from other organisms in the 
database. Moreover, we predicted the secondary struc- 
ture of this intron-H as shown in Figure 5. The second- 
ary structure modeling revealed distinctive IE sequences, 
namely, the non-expanding P5 and the expanding P9 
regions. 

Discussion 

To date, although a variety of introns from eukaryotes 
have been described in the rRNA gene loci of fungi [9], 
few introns in Phialophora species have been reported. 
An unusually small group 1 intron of 67 bps from the 
nuclear 18S rDNA has been described in a splicing 
study of Capronia semiimmersa, a teleomorph of P. 
americana which is known to be similar to P. verrucosa 
[20-22]. These small introns contain only PI, P7 and 
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Figure 4 A-C. - Diagrams for predicted secondary structure of P. verrucosa. [A]: intron-F from rDNA of PV1, [B]: intron-G from PV1 and [C]: 
intron-G from PV3. Capital letters indicate intron sequences and lowercase letters indicate flanking exon sequences. Arrows point to the 5' and 
3' splice sites. The guanosin cofactor-binding sites are marked with *. 
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lowercase letters indicate flanking exon sequences. Arrows point to 
the 5' and 3' splice sites. 



P10 elements, because most of the core regions com- 
mon in almost all other group 1 introns are missing. 
Four intron sequences have been reported or registered 
in dematiaceous fungi; namely, 283 bps within the small 
subunit (SSU) rDNA from Cadophora gregata f. sp. 
adzukicola [23], 339 bps within SSU from Cadophora 
finlandica (accession number: AF486119), 456 bps 
within the large subunit (LSU) rDNA from C. semiim- 
mersa [24] and 397 bps within LSU from Cladophialo- 
phora carrionii [24]. These introns have not been 
subjected to secondary structure analysis. Therefore, we 
aimed to identify the introns that we found in this study 
and to investigate the prevalence and phylogenetic rela- 
tionships of 28S group 1 intron at the intra-species level. 

The intron-F, G and H in the 28S rDNA of both spe- 
cies were found to belong to two subgroups, IC1 and IE, 
of group 1 intron. IC1 at L798 is the most common 
insertion position as shown in Table 1 and in the CRW 
website, and insertions at L1921 and L2563 were found 
comparatively in the database. The loss of most of P5 in 
the secondary structure of intron-H is believed to be a 
relatively recent evolutionary event [19]. The three 
insertions possessed all the ten elements (P1-P10) com- 
mon in group 1 introns. Enzymatic core regions are 
especially well conserved in primary and secondary 



structures, as described in previous reports [12,25], sug- 
gesting that they were derived from a common origin. 
Peripheral elements of the core have various forms and 
these variations have been used to subdivide introns 
into five major subgroups [17,26]. In this study, the phy- 
logeny obtained in Figure 2 and 3 showed that all IC1 
introns inserted into P. verrucosa have been surviving 
with base substitution/insertion/deletion, especially 
among peripheral elements as a consequence of some 
events after the individual insertion of IC1 at L798 and 
L1921, and may have spread by homing (e.g., [27-29]) or 
reverse splicing [30-32]. 

Comparisons of intron-F and G indicate comparative 
high sequence divergence within P. verrucosa wherein 
the sequence similarity among intron-F's was 94%, and 
99% among intron-G's with the exception of PV3 and 
90% among all the four intron-G's. The phylogenetic 
relationships among both introns were very similar to 
each other (Figure 2). These results suggest that they 
have diverged from a common origin. 

Thirty isolates of P. verrucosa and two isolates of P. 
americana possessed intron-F, G and H either individu- 
ally or as a combination of these introns. Genotypes 
based on the combinations of presence or absence of 
introns, type and position of insertions were established 
to discriminate among the isolates surveyed. As a result, 
five genotypes; namely, F, FG, FH, FGH and N were 
identified, as shown in Table 1. Type-F was isolated in 
all the countries where the strains used in the study. 
Intron distribution was found to have some correlation 
with geographic location, albeit the number of isolates 
used was small. For example, most of the Chinese iso- 
lates except for Yao-strain had only type-F. Isolates 
from South American continent had slight tendency to 
have an intron-H, wherein they were either type-FH or 
FGH. Intron-G's occurred as type-FG in the clinical iso- 
lates of Japan and China and as two type-FGH's in soil 
isolates of Brazil. In addition, according to an interpreta- 
tion from a different viewpoint, insight into possible 
correlation of geographic origin among introns from P. 
verrucosa strains have emerged from these insertion 
position results, namely, the spread of L798 among a 
large number of P. verrucosa isolates and the existence 
of L1921 and L2563 that coexist with the other intron 
insertions among the species and strains that have lost 
introns. L1921 positions are only seen in two clinical 
isolates from Japan and China and two isolates from 
Brazilian soil. The L2563 position seems to be specific 
to the South American continent in six environmental 
isolates. The possible correlation tendencies shown in 
our results have also been reported in previous studies 
described below. For example, group 1 intron CgSSU 
was found in the SSU rDNA of deuteromycetes mycor- 
rhizal fungus Cenococcum geophilum, and the intron- 
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positive isolates occurred mainly in North America and 
Europe and negative isolates in Western, Midwestern 
and Southern North America [33]. In other studies on 
intron distribution from a single species, four different 
group 1 intron combinations within LSU rDNA from 
entomopathogenic hyphomycete Beauveria bassiana 
were divided into 13 genotypes to investigate distribu- 
tion frequencies in the population and it was found that 
there was a tenuous correlation with geographic origin 
or insect host species [34]. Moreover, M. Marquez, et al. 
have found three intron insertion positions within LSU 
rDNA and established seven genotypes among 26 bio- 
control isolates for entomopathogenic anamorphic 
Metarhizum anisopliae [35]. 

Meanwhile, we found that five isolates of P. americana 
had no introns, even though two isolates were detected 
as type-F. Intron-loss strains might have been lost from 
taxa possessing intron-F in the common ancestor of the 
species during its evolution [36]. Further, the lateral 
transfers appear to have been rare events in P. ameri- 
cana. The fact that intron-F was found in almost all iso- 
lates of P. verrucosa, it is believed that intron-F may be 
specific to P. verrucosa. To confirm this hypothesis, 
more isolates are needed in the survey and the relation- 
ships of the clinical background of the individual 
patients and the ecological niches of saprobic isolates 
must be investigated. Further analysis of genotypes 
within the complete nuclear rDNA gene must be done 
and the presence of HE gene sequences must be ana- 
lyzed since they provide key information on intron phy- 
logeny and origin. This study is a first step in the study 
of introns in P. verrucosa and P. americana. 

Conclusion 

The three insertions within 28S rDNA of clinical and 
environmental isolates of P. verrucosa and P. americana 
allowed us to characterize them into five genotypes 
using agarose gel electrophoresis patterns. The two 
insertions, namely, intron-F and G, were characterized 
as subgroup IC1 by subjecting them to RT-PCR, sec- 
ondary structure and phylogenetic analysis to determine 
whether they are true introns, to characterize subgroup 
and to infer evolutionary relationships, respectively. 
Another insertion, intron-H, was characterized as an IE 
intron using BLAST search and by prediction of second- 
ary structure. Furthermore, we also developed a system 
to classify genotypes based on the presence and distribu- 
tion of group 1 introns and the distributions as DNA 
polymorphism among the two species. 

Methods 

Fungal strains and culture conditions 

We studied 34 P. verrucosa strains including of five clin- 
ical isolates as shown in Table 1. Seven P. americana 



strains including of three clinical isolates were used as 
allied species. All the isolates were preserved by using 
L-drying method and were sub-cultured on potato dex- 
trose ager (Difco) slant before extraction of genomic 
DNA. For an extraction of total RNA, liquid cultivation 
was performed in 50-ml Erlenmyer flask containing 20 
ml of potato dextrose medium at 30°C for seven days 
on a rotary shaker at 120 rpm. 

Extraction of genomic DNA and total RNA 

DNA extraction was performed using an InstaGene 
Matrix extraction kit (BioRad, Hercules, CA, USA) 
according to the manufacturer's instructions with minor 
revisions. Particularly, cells were ground with micro pes- 
tle before incubation at 56°C. The extracted DNA was 
then diluted 1:10 and used as template DNA for PCR 
amplification. 

Total RNA was extracted by using the Nucleic Acid 
Purification Kit MagExtractor (TM -RNA- TOYOBO, 
Osaka, Japan). The following procedures were done 
before carrying out the manufacturer's instructions. 
Approximately 20 mg (wet weight) of mycelia were 
washed with water and then rinsed with Schizosaccharo- 
myces pombe spheroplast buffer (20 mM citrate-phos- 
phate buffer (pH 5.6), 50 mM EDTA and 0.9 M 
sorbitol). This was followed by addition of 100 \i\ of buf- 
fer plus 20 units of Lyticase (L-5263; SIGMA, MO, 
USA) and 0.01 units of Chitinase (C-7809; SIGMA, MO, 
USA). The suspension was mixed by vortexing, incu- 
bated at 37°C for 30 min and added with 700 \i\ of a 
lysis-adsorption solution including 1% 2-mercaptoetha- 
nol and pricked with a pipette tip until the solution 
became viscous. The extracted RNA was treated with 
RNase-Free DNase Set (QIAGEN). Approximately more 
than 20 ng/|il RNA was obtained. 

PCR amplification and sequencing analysis 

A primer walking method was performed to obtain the 
sequences of the entire 28S rDNA region including ITS. 
PCR Master Mix (Promega, Madison, WI, USA) and 
TaKaRa LA Taq (TAKARA Bio Inc, Sigma, Japan) were 
used depending on the amplification sizes. The PCR 
conditions for PCR Master Mix consisted of denatura- 
tion for 4 min at 95°C, followed by 30 amplification 
cycles of denaturation at 94°C for 1 min, annealing at 
primer-dependent temperatures based on Tm values for 
1 min and extension at 72°C for 1.5 min, and then 1 
cycle of 5 min at 72°C. For TaKaRa LA Taq consisted of 
denaturation for 1 min at 94°C, followed by 30 cycles of 
denaturation at 98°C for 5 sec, annealing at primer- 
dependent temperatures for 30 sec and extension at 72° 
C for 2 min, and then 1 cycle of 72°C for 10 min. PCR 
products were purified with SAP-IT (USB Corporation, 
Cleveland, OH, USA) and then sequenced with primers 
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listed in Table 2 and the BigDye Terminator v3 Cycle 
Sequencing Kit (Applied Biosystems, Foster City, CA, 
USA) on an ABI Prism 3130 x 1 Sequencer (Applied 
Biosystems, Hitachi). The nucleotide sequences were 
determined from both strands. To determine base sub- 
stitutions and intron insertion positions, sequences were 
aligned by using the alignment function of GENETYX 
ver. 9.1.1 (GENETYX COOPERATION, Tokyo, Japan). 

Determining incidence of introns by agarose gel 

The extracted DNA was used as template DNA for the 
amplification of the insertion regions (intron-F, G and 
H). PCR was performed individually using PCR Master 
Mix and the primer pair inF-F and inF-R for intron-F 
and inG-F and inG-R for intron-G which we newly 
designed. Primer pair L2563F and L2563R for intron-H 
was designed based on sequences of exon and group 1 
intron on CRW website, because the intron was not 
inserted in the five representative strains used. PCR con- 
ditions were the same as described above and the result- 
ing DNA fragments were resolved by electrophoresis on 
a 2% agarose gel (NuSieve® 3:1 Agarose, TAKARA Bio 
Inc, Sigma, Japan) in Tris-borate-EDTA buffer. Presence 
or absence of individual intron was listed as positive/ 
negative in Table 1. In addition, the strains were cate- 
gorized into five intron types; namely, F, FG, FH, FGH 
and N on the basis of the intron insertions. 

RT-PCR and colony sequencing 

The RT-PCR from total RNA was performed using a 
Superscript ™ III One Step RT-PCR System with Plati- 
num Taq DNA Polymerase (Invitrogen, CA, USA) 
according to the manufacturer's instructions. Ten 
microliters of RNA was used for a 50 \A reaction mix- 
ture consisting of 25 \A of ReactionMix, 2 \A of Super- 
Script™ III RT/Platinum Taq Mix, 1 ul of each primer 
inF-F (10 uM) and inF-R (10 uM) for RT-PCR of 
intron-F and 11 \i\ of autoclaved distilled water. For RT- 
PCR of intron-G, primer pair inG-F and inG-R was 
used. RT-PCR was carried out in the following condi- 
tions: cDNA synthesis at 55°C for 30 min, denaturation 
at 94°C for 2 min, and PCR amplification at 40 cycles of 
94°C for 15 sec, 55°C for 30 sec and 68°C for 1.5 min 
and final extension at 68°C for 5 min. Amplification 
products were eluted in 3.5% polyacrylamide gel in tris- 
acetate-EDTA buffer on an electrophoresis run condi- 
tion of 100 V for 30 min and followed by 75 V for 25 
min, together with genomic DNAs amplified with the 
same primer pairs as control (shown in Figure 1). 

The RT-PCR products were purified with the SUPREC- 
PCR (TAKARA Bio Inc, Sigma, Japan) and ligated into 
the pGEM-T Easy Vector System (Promega, Madison, 
WI, USA). Plasmids were transformed into E. coli compe- 
tent cells (ECOS TM Competent E. coli, JM109, NIPPON 



GENE Co., LTD., Japan). Transconjugants were selected 
on LB agar plates containing 50 (ig/ml ampicilin and 40 
(ig/ml of 5-bromo-4-chloro-3-indoyl-P-D-galactopyrano- 
side (X-Gal). The presence of the expected insert was 
confirmed by PCR and agarose gel electrophoresis. The 
inserts were sequenced with T7 (5'-TAATACGACT- 
CACTATAGGG-3') and M13 reverse primers (5'- 
AGGAAACAGCTATGACCATGA-3'). 

Phylogenetic analysis of introns from P. verrucosa 

Nucleotide sequences were aligned using the BioEdit 
program version 7.0.9.0 [37]. For phylogenetic analysis, 
alignment gaps were treated as missing data and ambig- 
uous positions were excluded from the analysis. NJ ana- 
lysis [38] as distance matrix method and MP analysis as 
character state method were carried out using PAUP 
4.0bl0 [39]. For NJ analysis, the distances between 
sequences were calculated using Kimura's two-para- 
meter model [40]. MP analysis was undertaken with the 
heuristic search option using the tree-bisection-recon- 
struction (TBR) algorithm with 1000 random sequence 
additions to find the global optimum tree. All positions 
were treated as unordered and unweighted. The maxi- 
mum tree number was set at 10 4 . To estimate clade 
support, the bootstrap procedure of Felsenstein [41] was 
employed with 1000 replicates in both MP and NJ ana- 
lyses. Bootstrap (BS) values higher than 50% are 
indicated. 

Alignment and phylogenetic analysis of core sequences 

For the comparison with highly conserved sequences of 
subgroup IC1 from 20 taxa, sequences of elements of P, 
Q, R and S and the pairing segment P3 were obtained 
from DDBJ database (accession numbers shown after 
sample name in Figure 3). These regions do not include 
IGS, because the sequences in the upstream region of 
intron insertion positions do not share a common IGS 
[42]. The NJ tree was constructed after alignment of all 
the sequences, which ranged from 57 to 60 bps (Addi- 
tional File 2). Insertion positions are shown after the 
sample ID or accession number. The insertion position 
numbering of the taxa refers to the 23S nucleotide 
sequence of E. coli [10] except for the Exophiala cali- 
cioides which is based on the 18S of Saccharomyces cer- 
evisiae. Taxonomic affiliation was indicated by letters in 
parentheses; namely, [A], Fungi/Ascomycota; [Ac], 
Acanthamoebidae; [C], Chlorophyta; [H], Heterolobosea; 
[M], Mycetozoa and [R], Rhodophyta. 

Secondary structure modeling 

The secondary structures are proposed from modeling by 
Michel et al. [14,26,43] and computational analysis was 
done using the Mfold web server available at http:// 
mfold.rna.albany.edu/ [44] and GENETYX Ver.9 software, 
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with manual adjustments. The pairing segments of Pl- 
P10 locations are indicated in Figure 4 and 5. Moreover, 
the model was manually optimized based on previous 
studies of group 1 introns [17,45-47]. 

Additional material 
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Additional file 1: Schematic representation of the large ribosomal 
subunit 28S gene. The hatched and dotted boxes correspond to the 
group 1 intron of P. verrucosa inserted at positions 798, 1921 and 2563 
relative to the 23S rDNA of the £ coli J01965 sequence. The numbering 
in the parentheses is relative to the ITS and 28S rDNA sequence of P. 
verrucosa. 

Additional file 2: Partial alignment of IC1 introns of P. verrucosa 
and selected introns from the database. Highly conserved sequences 
of the elements of P, Q, R and S and the pairing segment P3 are also 
shown. Intron insertion positions relative to £ coli are given after the 
sample ID or taxon name. * indicates the insertion position relative to 
the 18S rDNA of the 5. cerevisiae sequence. Letters in parentheses 
indicate taxonomic affiliation: [A], Fungi/Ascomycota; [Ac], 
Acanthamoebidae; [C], Chlorophyta; [H], Heterolobosea; [M], Mycetozoa; 
[R], Rhodophyta. 

Additional file 3: Alignment of intron-F used for the phylogenetic 
analysis and the modeling of secondary structure. The gaps were 
marked with dashes. The highly conserved (ribozymatic core) regions of 
the P, Q, R and S were marked with dotted lines. Boxed nucleotides 
participate in the pairing segments of PI -PI 0 of the secondary structure 
model. 

Additional file 4: Alignment of intron-G used for the phylogenetic 
analysis and the modeling of secondary structure. The gaps were 
marked with dashes. The highly conserved (ribozymatic core) regions of 
the P, Q, R and S were marked with dotted lines. Boxed nucleotides 
participate in the pairing segments of PI -PI 0 of the secondary structure 
model. 



Takizawa et al. BMC Microbiology 201 1, 1 1:94 Page 1 2 of 1 2 

http://www.biomedcentral.eom/1 471 -21 80/1 1 /94 



30. 



31. 



32. 



33. 



34. 



35. 



36. 



37. 



39. 



40. 



41. 



42. 



43. 



44. 



45. 



46. 



47. 



49. 



50. 



Woodson SA, Cech TR: Reverse self-splicing of the Tetrahymena group 1 
intron: Implication for the directionality of splicing and for intron 
transposition. Cell 1989, 57(2):335-345. 

Roman J, Woodson SA: Reverse splicing of the Tetrahymena IVS: evidence 
for multiple reaction sites in the 23S rRNA. RNA 1995, 1:478-490. 
Roman J, Woodson SA: Integration of the Tetrahymena group 1 intron 
into bacterial rRNA by reverse splicing in vivo. Proc Natl Acad Sci USA 
1998, 95:2134-2139. 

Shinohara ML, LoBuglio KF, Rogers SO: Group-1 intron family in the 
nuclear ribosomal RNA small subunit genes of Cenococcum geophilum 
isolates. Curr Genet 1996, 29(4):377-387. 

Wang C, Li Z, Typas MA, Butt TM: Nuclear large subunit rDNA group 1 
intron distribution in a population of Beauveria bassiana strains: 
phylogenetic implications. Mycol Res 2003, 1 07(1 0):1 189-1200. 
Marquez M, Iturriaga EA, Quesada-Moraga E, Santiago-Alvarez C, Monte E, 
Hermosa R: Detection of potentially valuable polymorphisms in four 
group 1 intron insertion sites at the 3'-end of the LSU rDNA genes in 
biocontrol isolates of Metarhizium anisopliae. BMC Microbiol 2006, 6:77-84. 
Haugen P, Simon DM, Bhattacharya D: The natural history of group I 
introns. Trends Genet 2005, 21 (2):1 11-119. 

Hall TA: BioEdit: a user-friendly biological sequence alignment editor and 
analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 1999, 
41:95-98. 

Saitou N, Nei M: The neighbor-joining method: a new method for 
reconstructing phylogenetic trees. Mol Biol Evol 1 987, 4:406-425. 
Swofford DL: PAUP: Phylogenetic analysis using parsimony [and other 
methods]. Sinauer, Sunderland, MA; 2003, Version 4. 
Kimura M: A simple method for estimating evolutionary rates of base 
substitutions through comparative studies of nucleotide sequences. J 
Mol Evol 1980, 16(2):111-120. 

Felsenstein J: Confidence limits on phylogenies: An approach using the 
bootstrap. Evolution 1985, 39(4)783-791. 

Hoshina R, Imamura N: Phylogenetically close group 1 introns with 
different positions among Paramecium bursaria photobionts imply a 
primitive stage of intron diversification. Mol Biol Evol 2009, 

26(6):1309-1319. 

Cech TR, Damberger SH, Gutell RR: Representation of the secondary and 
tertiary structure of group 1 introns. Nat Struct Biol 1994, 1 (5):273-280. 
Zuker M: Mfold web server for nucleic acid folding and hybridization 
prediction. Nucleic Acids Res 2003, 31 (1 3):3406-341 5. 
Johansen S, Johansen T, Haugli F: Structure and evolution of myxomycete 
nuclear group 1 introns: a model for horizontal transfer by intron 
homing. Curr Genet 1992, 22:297-304. 

Egger K: Sequence and putative secondary structure of group 1 introns 
in the nuclear-encoded ribosomal RNA genes of the fungus 
Hymenoscyphus ericae. Biochimica et Biophysica Acta (BBA) - Gene Structure 
and Expression 1995, 1261(2):275-278. 

Perotto S, Nepote-Fus P, Saletta L, Bandi C, Young JPW: A diverse 
population of introns in the nuclear ribosomal genes of ericoid 
mycorrhizal fungi includes elements with sequence similarity to 
endonuclease-coding genes. Mol Biol Evol 2000, 17(1):44-59. 
White TJ, Bruns TD, Lee SB, Taylor JW: Amplification and direct 
sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR 
Protocols: a Guide to Methods and Applications. Edited by: Innis MA, Gelfand 
DH, Sninsky JJ, White TJ. London: Academic Press; 1990:315-322. 
O'Donnell K: Fusarium and its near relatives. In The fungal holomorph: 
mitotic, meiotic and pleomorphic speciation in fungal systenatics Edited by: 
Reynolds DR, Taylor JW 1993, 225-233. 

McCullough MJ, Clemons KV, Stevens DA: Molecular and phenotypic 
characterization of genotypic Candida albicans subgroups and 
comparison with Candida dubliniensis and Candida stellatoidea. J Clin 
Microbiol 1999, 37(2):41 7-421. 

Williamson ECM, Speers D, Arthur IH, Harnett G, Ryan G, Inglis TJJ: 
Molecular epidemiology of Scedosporium apiospermum infection 
determined by PCR amplification of ribosomal intergenic spacer 
sequences in patients with chronic lung disease. J Clin Microbiol 2001, 
39(1):47-50. 



doi:1 0.1 1 86/1 471 -21 80-1 1 -94 

Cite this article as: Takizawa et al:. Occurrence and characteristics of 
group 1 introns found at three different positions within the 28S 
ribosomal RNA gene of the dematiaceous Phialophora verrucosa: 
phylogenetic and secondary structural implications. BMC Microbiology 

2011 11:94. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



